Anomalous Temperature Dependence of Lower Critical Field in Ultraclean URu2Si2 
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To investigate a mysterious superconducting state of URu2Si2 embedded in the so-called hidden 
order state, the lower critical field Hci is precisely determined down to 55 mK for Ji"|| a and H\\ 
c. For this purpose, the positional dependence of the local magnetic induction is measured on 
ultraclean single crystals (Tc = 1.4 K) with residual resistivity ratio exceeding 700. We find that 
the temperature dependence of Hci significantly differs from that of any other superconductors. 
The whole HdiT) for H\\ a are well explained by the two superconducting gap structures with line 
and point nodes, which have been suggested by the recent thermal conductivity and specific heat 
measurements. On the other hand, for H\\ c, a change of slope with a distinct kink in Hci{T), which 
cannot be accounted for by two gaps, is observed. This behavior for H\\ c sharply contrasts with 
the cusp behavior of Hci (T) associated with a transition into another superconducting phase found 
in UPts and Ui-ajThj^Beia. The observed anomalous low-field diamagnetic response is possibly 
related to a peculiar vortex dynamics associated with chiral domains due to the multicomponent 
superconducting order parameter with broken time reversal symmetry. 



I. INTRODUCTION 

Exotic superconducting state with nontrivial Cooper 
pairing in heavy fermion systems continues to be a cen- 
tral focus of investigations in strongly correlated elec- 
tron systems.— 1^ Among them, URu2Si2 has mystified 
researchers since the heavy fermion superconductivity 
occurs deep inside the mysterious 'hidden order' state 
whose transition temperature is Th = 17.5 Kj^"— 

Several salient features of the superconducting state 
have been reported in URu2Si2. According to several 
experimental observations, most of the carriers (~ 90 %) 
disappear below Thr^^ resulting in a semimetallic elec- 
tronic structure with a density one order of magnitude 
smaller than in other heavy fermion superconductors. 
Superconductivity with such a low density is remark- 
able since the superfluid density is very low in some way 
reminiscent of underdoped cuprates; the superconduc- 
tivity by itself is an exceptional case of pairing among 
heavy electrons with a long Fermi- wavelength in a nearly 
semimetallic system. Moreover, pressure studies revealed 
that the superconductivity coexists with the hidden or- 
der phase having no intrinsic magnetic order but cannot 
coexist with the antiferromagnetic order with large mag- 
netic moment, in contrast with other heavy fermion 
compounds where the superconductivity often coexists 
with magnetic orders. Thus, although the genuine hid- 
den order parameter is still an open question, the hidden 
order may provide an intriguing stage for a new type of 
unconventional superconducting state. 

Despite these studies, very little is known about the 
exotic superconducting state of URu2Si2, mainly be- 
cause the superconductivity is extremely sensitive to 



disorder .A3. Recent studies using ultraclean single crys- 
tals with very large residual-resistivity-ratio (RRR) 
have enabled us to develop an advanced understand- 
ing about the superconducting state of URu2Si2i^'^ 
A number of unprecedented superconducting properties 
have been reported, including superconductor-insulator- 
like first-order transition at upper critical field Hc2^ 
flux line lattice melting at sub-Kelvin temperatures^^ 
and quantum transport of quasiparticlesJ^ Moreover, 
it has been shown that below Th the nearly perfect 
compensation (equal number of electrons and holes) 
is realized.— According to the quantum oscillation,—""— 
electronic transport?^ specific healii^ and thermal con- 
ductivity measurements?^ the presence of the light spher- 
ical hole band and anisotropic heavy electron band has 
been suggested. This indicates that URu2Si2 is an essen- 
tially multiband superconductor. Recent angle-resolved 
thermal conductivity and specific heat measurements 
have suggested two distinct superconducting gap struc- 
tures having different nodal topology with horizontal line 
nodes in the hole band and point nodes in the electron 
hand^'^— From the group symmetry analysis, a chiral 
d-wave state with a form. 



A(fe) = Aak^ikj: ± iky 



(1) 



has been proposed, which breaks time reversal symmetry 
(TRS)^'ia^ 

Up to some years ago, the only known condensate 
with broken TRS was the A-phase of superfluid ^He. 
More recently, several candidates of unconventional su- 
perconducting phases which break TRS have been re- 
ported in several classes of strongly correlated materi- 
als, including ruthenate Sr2Ru04,^^ filled skutterudite 
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PrOs4Sbi2r^ and U-based heavy fermion compounds, 
UPtg^ and Ui-j^Thj-Beisi^^ In these materials, muon 
spin relaxation measurements in zero field reveal the 
development of spontaneous magnetic moment below 
Tc and in Sr2Ru04, nonzero Kerr rotation has been 
reportedi^ Theoretically, it is predicted that in a super- 
conductor with multicomponent order parameter, chiral 
domain walls separating different degenerate supercon- 
ducting states are formed i^"— Then the system should 
spontaneously generate supercurrents at the edge and the 
domain walls. Moreover, it is proposed that these chiral 
domains should strongly influence the vortex penetra- 
tion especially at low fieldS)^ These chiral domains and 
low-field magnetic response in unconventional supercon- 
ductors with broken TRS have been explored by several 
experimental techniques, '^^^SS, however there is no direct 
evidence of such anomalies in the supercurrents as well 
as vortex penetration. 

In this paper, to shed further light on the exotic su- 
perconducting properties of URu2Si2, we investigated the 
low-field diamagnetic response of ultraclean single crys- 
tals. The results provide strong evidence for the highly 
unusual superconducting state of URu2Si2 embedded in 
the hidden order state. 



II. EXPERIMENTAL 

High-quality single crystals of URu2Si2 with RRR = 
700 were grown by the Czochralski pulling method in 
a tetra-arc furnace.— The well defined superconducting 
transition was confirmed by the specific heat measure- 
ments. Experiments have been performed on single crys- 
tals with typical dimensions of 2.3 x 0.75 x 0.15 mm'^ (see 
the inset of Fig. 3) down to 55 mK by using a dilution 
refrigerator. 

The standard technique to determine the lower criti- 
cal field Hci is to measure the dc magnetization in the 
bulk crystals. However, reliable determination of Hd is 
a difficult task in the presence of the flux pinning and 
magnetic relaxation. To avoid these difficulties, we de- 
termined Hci by directly detecting the positional depen- 
dence of the field Hp at which flux penetration occurs 
from the edge of the crystal.— To measure the local mag- 
netic induction, we used micro Hall sensors tailored in a 
GaAs/AlGaAs heterostructurei^ 

First we studied the pinning properties in the crystals 
by measuring the magnetic field distribution in the vortex 
state, which is obtained by the scanning micro Hall sen- 
sor with an active area of 2 x 2 jixn^ . The Hall sensor was 
scanned on the surface area of the crystal by using the 
piezoelectric device. The distance between the Hall sen- 
sor and the crystal surface was kept constant by monitor- 
ing the tunneling currents. Next we measured the local 
magnetic induction very precisely by using a miniature 
Hall-sensor array with each active area of 5 x 5 ^jp? 
(the center-to-center distance of neighboring sensors is 
20 /im). The crystal was directly placed on top of the 
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FIG. 1: (color online). The positional dependence of the local 
magnetic induction B near the edge region of the crystal, 
which is determined by the scanning Hall-probe microscopy. 
The dashed line is the position of the edge (For details, see the 
text.) (a) The profile of magnetic induction in the Meissner 
state. The external magnetic field {fioH = 0.05 mT) is applied 
parallel to the c axis, (b) The profile of magnetic induction 
in the vortex state at several fields at T = 0.5 K. (c) Contour 
image of (b). 



array (see the inset of Fig. 2(a)). 

In both measurements, earth magnetic field was 
shielded by /i-metal. The residual field at the sample 
space is less than 0.5 /zT. In all measurements, the ex- 
ternal field was applied after the sample was zero-field 
cooled to the desired temperature from the temperature 
above Tc. 



III. RESULTS 

Figure 1 (a) displays the profile of the magnetic induc- 
tion near the edge of ab plane in the Meissner state when 
the external magnetic field ^j-qH = 0.05 mT is applied 
parallel to the c axis. The profile is measured by the 
scanning Hall-probe microscopy. With approaching from 
the outside of the crystal, the magnetic induction is en- 
hanced to B ^ 0.08 mT followed by a sharp drop near the 
edge of the crystal. The induction is perfectly screened 
inside the crystal. Here we defined the edge at the peak 
of the local magnetic induction, as shown by the dotted 
line. The enhancement of the magnetic induction at the 
edge is due to the field induced by the Meissner shielding 
currents fiowing the edge region^SS, The perfect screening 
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FIG. 2: (color online), (a) Local magnetic induction Bodge 
measured by the sensor at the edge of the crystal (Sensor 2) , 
as a function of H for a. The triangles indicate the flux 
penetration field Hp. Each data is vertically shifted for clarity. 
Inset: Photograph of the miniature Hall-sensor array, and a 
schematic illustration of the experimental setup. Sensor 1 
locates just outside of the crystal, (b) The same plot for 
H\\ c. Inset is the positional dependence of Hp. 
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FIG. 3: (color online). Inset: Photograph of the single crystal 
URu2Si2 used in the present study. Local magnetic induction 
is measured at the edge regions A, B, and C. Main panel: 
Temperature dependence of the lower critical fields Hci for 
H\\ a (solid squares) and H\\ c (solid circles), which are de- 
termined at the edge region B shown in the inset. Arrows 
indicate a kink anomaly at Tq = 1.2 K and a cusp behav- 
ior at Tp ~ 1 K. The dotted lines are the fits to the linear 
dependence. Open circles and open squares show H^i{T) de- 
termined at the edge regions A and C shown in the inset, 
respectively. 



is achieved in a short distance at d > 14 /im just below Tc, 
where d is the distance from the edge. This demonstrates 
the homogeneous superconducting state of the crystal. 

Figures 1(b) and (c) display the distribution of the lo- 
cal magnetic induction at several external magnetic fields 
associated with the flux penetration at T = 0.5 K. The 
magnetic induction decays monotonically with the dis- 
tance d from the edge. This indicates the Bean critical 
state dominated by the bulk flux pinning. In this case, 
Hci is determined from the local magnetic induction near 
the edge, at which the first flux penetration is most sensi- 
tively detectedi^ This is in contrast to the case when the 
geometrical surface barrier (determined by sample shape) 
is important rather than bulk pinning, in which the mag- 
netic flux distribution shows a dome-like shape and the 
flux penetration field can be detected at the center of the 
crystal. 

The positional dependence of the local magnetic induc- 
tion i?odgc were measured accurately by using the Hall- 
sensor array near the edge for H\\ a and H\\ c. The 
sample edge is located between the Hall sensors 1 and 2 
(see the inset of Fig. 2), which is confirmed by the Meiss- 
ner response of each sensor, as shown in Fig. 1(a). In 
fact, sensor 2 exhibits a perfect Meissner response, while 
slight enhancement of B from fioH is observed in sensor 
1. Figures 2(a) and (b) display the field dependence of 
Badge measured by sensor 2. The flux penetration fields 
Hp shown by the triangles are clearly resolved by the de- 



viation from the Meissner state (Sedge — 0). The inset of 
Fig. 2(b) depicts the positional dependence of Hp. Obvi- 
ously Hp is position-independent at sensors 2 and 3 which 
locate close to the edge (d < 40 /xm) even at low temper- 
atures. On the other hand. Hp is enhanced at sensors 4 
and 5, which locate at d > 40 //m, at low temperatures. 
These results again represent the Bean critical state in 
this crystal and assure that we can accurately determine 
Hci from Hp measured by sensors near the edge. 

Figure 3 depicts the temperature dependence of H^^ 
and H^i , where H^i and H^i are the lower critical fields 
for a {H\\ c), respectively. We evaluate H^i = 
3.82Hp and H^i — 1.15Hp by using the expression 



tanh ■\/0.366/a' 

where a and b are the width and the thickness of the 
crystal, respectively, for a platelet sample. This ex- 
pression takes into account the geometrical effect in an 
nonellipsoidal shape, where the demagnetization effect 
in ideal ellipsoids is modified.— The magnetic penetra- 
tion depth is evaluated from Hci to Aa — 0.8 /im, 
which is quantitatively consistent with the /xSR results of 
Aa = 0.7—1 /im.— According to the Ginzburg-Landau 
theory, the temperature dependence of Hd anisotropy, 
•jH^i = Hci/ Hci, t>e different from that of Hc2 in 
multiband superconductors.^'^''*^ At T Tc, however. 



4 



7_ff^i should coincide with the Hc2 anisotropy. The ex- 
perimentahy determined JHci — 3 close to is then 
fully consistent with the reported Hc2 anisotropy near 
Tcr^ giving us confidence in the accuracy of the result. 
As seen in Fig. 3, H^i increases linearly below ^ 0.8 K 
with decreasing temperature. A cusp structure can be 
seen at Tp =1 K. Meanwhile, H^^ increases steeply be- 
low Tc and increases linearly after exhibiting a distinct 
change in the slope with a kink at Tq = 1.2 K. As shown 
by the dotted lines, H^^ at low temperatures increases 
linearly down to 55 mK, while H^^ exhibits a tendency 
to saturation below 200 mK. 



IV. DISCUSSION 

First we discuss the T-dependence of Hd in the low 
temperature regime. Since Hd is proportional to the su- 
perfiuid density, the T-linear dependence of down 
to low temperatures indicates the presence of line nodes 
in the superconducting gap function. Moreover, the ten- 
dency to saturation in H^^{T) below 200 mK indicates 
that the line nodes are located parallel to the ah planes 
(horizontal node). This is because the supercurrents al- 
ways flow parallel to the nodal planes for H\\ c, while for 
1? II a supercurrents have a component which flows across 
the horizontal nodes. This component is perpendicular 
to the velocities of the quasiparticles around the nodes, 
which reduces the contributions of nodal quasiparticles 
to the superfiuid density. Since in the multiband super- 
conductors the penetration depth and hence the lower 
critical field are governed by the band with light mass, 
it is natural to consider that the horizontal line node 
locates in the spherical light hole band. This is consis- 
tent with the previous results of angle-resolved thermal 
conductivity3^i2iL 

Most remarkable features are anomalies at Tq for H\\ c 
and Tp for a. In the specific heat measurements, no 
anomaly has been observed at Tp in the present crystals 
with large RRR values. However, a peak structure in 
the heat capacity has been reported in the vicinity of 
Tp for crystals with RRR < 100, possibly due to the 
inhomogeneous distribution of Tc<^ Therefore we cannot 
rule out a possibility that the observed Tp-anomaly in H^^ 
is due to a tiny portion with low-Tc phase in the crystal. 
On the other hand, no anomaly has been reported in 
the heat capacity at Tq for any crystals, indicating no 
evidence of the low-Tj, phase as an origin for the Tq- 
anomaly. It is also highly unlikely that the anomaly at 
Tq arises from the inhomogeneous flux penetration of 
the crystal, because of the following reasons. First, the 
open squares and circles in Fig. 3 show H^^ determined at 
different edge regions in the same crystal. The distinct 
kink anomaly at Tq is quite well reproduced. Second, 
as shown in Fig. 1(a), the Meissner shielding currents 
flow homogeneously and the crystal has well deflned edge. 
These results lead us to conclude that the anomaly of 
Hd{T) at Tq is intrinsic. We stress that the anomaly 
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FIG. 4; (color online). Comparison of h""'^ = 
^^S"(r)/i?S"(0) (symbols) with (a) ^/TiJTif and (b) nj, re- 
spectively. The dotted and dashed lines indicate the super- 
fluid densities for the hole band with Aj, = LSksTc and the 
electron band with Ae — A.OksTc, respectively, which are ob- 
tained by assuming the illustrated nodal topologies suggested 
by the thermal conductivity measurementsi^ The solid thick 
and thin lines represent the results of the two-gap fltting and 
the calculations for other gap symmetries, respectively. 



at Tq has never been reported in the crystals with lower 
RRR values4^ 

Now we try to fit the temperature dependence of 
H^^ and H^^ in accordance with the multiband model. 
Recent angle-resolved thermal conductivity and specific 
heat measurements revealed two distinct superconduct- 
ing gap structures having different nodal topology with 
horizontal line nodes in the hole band and point nodes in 
the electroni^ d^i^° To see whether multiband effect can 
explain anomalous behaviors of the lower critical fields, 
we discuss the temperature dependence of and 
in terms of the two-gap model below. In the two-band 
model, the in-plane and out-of-plane superfiuid density 
normalized by their values at T — OK, n° and n^, re- 
spectively, can be written as 

nliT) = x^nUT) + {l-x^)nl{T), i = a,c (3) 
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where nj^ and n\ are the normahzed superfluid density of 
hole and electron bands, respectively, is the ratio of 
the electron and hole contributions given by 



v-i _ 



kF,h(e)l 



(4) 
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where Vp ^^^^^ is the Fermi velocity of holes (electrons^ 
The hole mass determined by the dHvA measurements is 
nearly isotropic « « ISmor^*^ where {nT-h) is 
the mass of the hole along the a (c) axis and mg is the free 
electron mass. The heavy electron mass is anisotropic 
and is estimated to be « SStoq and « SOStoq 
from the large Sommerfeld coefficient in the heat capac- 
ity (-y ~ 80mJ/K^mol)^ and the anisotropy of upper 
critical field Hc2-'^ Then we obtain x"' « 0.87 and x'^ 
0.95. Since theses values are close to unity, the superfluid 
density, particularly its temperature dependence at low 
temperatures, is governed by the hole band. This again 
provides a strong support to the horizontal line nodes lo- 
cated on the hole bands as mentioned before. Moreover, 
the presence of point nodes has been suggested along the 
c axis in the heavy electron bandsi^ii^ Then, we calcu- 
late the superfluid density by assuming the gap functions 
Ah{k,T) — Ah{T) X 2 sin 6* cos 6' with point and horizon- 
tal line nodes for hole bands and Ag {k,T) — Ag (T) x sin 9 
with c-axis point nodes for electron bands. 

The normalized lower critical flelds, /i°'^(T) = 
H'^i {T) / H'^i {Q) , are related to the superfluid density 
as h" = n1 and h°- = ^nfrij. First we try to fit h°-{T). 
Because the hole band dominates at low temperatures, 
Ah{k,T) can be determined by /i"(T). The best fit is 
obtained by A^(0) = IMbTc and Ae(0) = iMsTc. As 
shown by the thick line in Fig. 4(a), the fitting result 
reproduces well the overall temperature dependence of 
/i°(r), including the observed tendency toward satura- 
tion at low temperatures. We note that this Ah{0) value 
is close to A/i(0) ^ hvF,h/'n'£.h ^ l.SfcsTc, which is ob- 
tained from vp^h ^ 2 X 10'*m/si^ and ^ 25 nm. Here 
~ ^ ^o/2Trfj,QH^2 is the coherence length of the hole 
band estimated from the "virtual upper critical field" of 
the hole band A*oi^c2 ^ 0.5 T.^ In Fig. 4(b), h%T) cal- 
culated by using the same A/i(0) and Ae(0) values is 
plotted by the thick line. In sharp contrast to h°'{T), the 
calculation strongly deviates from the data. We tried to 
fit the data by assuming various A;i(0) and Ae(0) values 
and other gap symmetries, but could not reproduce the 
data, particularly the anomaly at Tq. 

Based on these results, we conclude that the multigap 
effect cannot be an origin of the anomaly at Tq. Al- 
though a change in the slope of Hci{T) has been reported 
in UPta'''^ and (Ui_:rTh^)Bei3^i^ there are crucial dif- 
ferences from our observation. First, in URu2Si2 Hci{T) 
is strongly suppressed below Tq, while in UPts and 
(Ui_a:Th2:)Bei3 Hci{T) is enhanced below the kink tem- 
perature. Second, in URu2Si2 the kink anomaly is ob- 
served solely for H\\ c, while in UPta and (Ui_a;Tha;)Bei3 
the anomaly is observed in any field directions. In the lat- 



ter compounds, the kink anomaly has been attributed to 
the transition into another superconducting phase with 
different gap structure. 

The observed strong suppression of H^i{T) indicates 
the drastic change of the vortex entry to the sample at 
Tq. Here we try to explain qualitatively this anoma- 
lous behavior of Hci{T) in terms of a peculiar vortex 
dynamics associated with chiral domains due to the niul- 
ticoniponent superconducting order parameter with bro- 
ken TRS. It is natural to consider that the superconduct- 
ing order parameter couples to the hidden order, because 
the superconductivity appears only in the hidden order 
phase as shown by recent pressure studies>ii To discuss 
the pairing state in the presence of the hidden order, we 
investigate the dimensionless Ginzburg-Landau free en- 
ergy which describes the pairing state with kz{kx + iky) 
form at zero magnetic field and low temperature. 



. T-T, 



Te T, 



+ /33imrh2r 

Here the order parameter is written as 

A,,,(fc)=ryi0-\(fe)-Hr;2C,,(fe), 



(5) 



(6) 



where (j)^\{k) and (f^\{k) are the order parameters hav- 
ing kzkx and k^ky symmetry, respectively. When the sys- 
tem has a four-fold rotation symmetry in the ab plane, 
these pairing states must have the same transition tem- 
perature, and Tp' = Tc. Then, the TRS is spontaneously 
broken just below T^ within the BCS theory since /?2 > 
and — 2/32 < 0. One the other hand, the four-fold 
rotation symmetry may be broken by the hidden order 
which occurs at Th = 17.5 K. When this is true, the tran- 
sition temperature of the kzkx state differs from that of 
the kzky state, namely T^ ^ Tc. (We can assume T^. < Tc 
without loss of the generality.) Then, the pairing state 
without broken TRS (ryi, 772) oc (1, 0) is stable just below 
Tc and the second-order phase transition to the pairing 
state (771,772) oc (Ijiia) occurs at 



Tc2 — Tc — 



Tc - n 

2/32 - fiz 



<Tc. 



(7) 



The TRS is broken at this second superconducting tran- 
sition. 

Assuming Tc2 — Tq, the order parameter with broken 
TRS can produce domain walls along the c axis at Tq. 
As suggested in Ref. 30, magnetic fields then can pene- 
trate inside through the domain wall even in the Meissner 
state. Thus the formation of domain walls gives rise to 
a peculiar reduction of the penetration field only for H 
II c. As shown in Figs. 1(a), (b) and (c), our scanning 
Hall probe detects no apparent inhomogeneous field dis- 
tribution, which implies that such domains, if exist, are 
smaller than the size of Hall sensor 2/xm). We also 
note that the kink anomaly at Tq has never been ob- 
served in the crystals with low RRR values?^ implying 
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that the domain formation is sensitive to the impurities. 
It should be noted that a possible presence of such a 
domain structure, which influences the vortex dynamics, 
has been suggested in UPts, Ui_2;Tha;Bei3, Sr2Ru04 and 
PrOs4Sbi2^^ 

The second superconducting transition discussed here 
is similar to the transition from the p-wave pairing state 
with TRS to the chiral p-wave state without TRS which 
has been proposed for Sr2Ru0442i The four-fold rota- 
tion symmetry is broken by the in-plane magnetic field 
in Sr2Ru04, while the hidden order may spontaneously 
break the four- fold rotation symmetry in URu2Si2 with- 
out applying any external field. We note that local four- 
fold symmetry breaking in the hidden order phase has 
been recently proposed theoretically— 

V. CONCLUSION 

We have determined accurately the lower critical fields 
of ultraclean URu2Si2 by the positional dependence of 
the local magnetic induction. The lower critical fields 
exhibit several distinct features, which have never been 
observed in any other superconductors. The tempera- 



ture dependence of Hd at low temperatures suggests the 
horizontal line nodes in the energy gap in the light hole 
band. We show that the whole Hci{T) ioi H\\ a can 
be explained by the two gaps with line and point nodes, 
which is consistent with the previous reports. In sharp 
contrast, for H\\ c we find a distinct kink in the slope of 
Hci{T) at 1.2 K, which cannot be accounted for by the 
two gaps. This anomalous low-field diamagnetic response 
has been discussed in the light of the multicomponent 
order parameter with broken TRS. The highly unusual 
superconducting state of URu2Si2 deserves further theo- 
retical and experimental investigations. 
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